Layered Structures Favor Superconductivity in Compressed Solid SiH4 
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The electronic and lattice dynamical properties of compressed solid SiH4 have been calculated in 
the pressure range up to 300 GPa with density functional theory. We find that structures having a 
layered network with eight-fold SiHs coordination favor metallization and superconductivity. SiH4 
in these layered structures is predicted to have superconducting transition temperatures ranging 
from 20 to 80 K, thus presenting new possibilities for exploring high temperature superconductivity 
in this hydrogen-rich system. 
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SiEU offers a source for high purity silicon in epitaxial 
and thin film deposition which is at the base of electronics 
and microdevices. There is ongoing interest in this ma- 
terial as well due to the suggestion of Ashcroft [l| that 
SiEU would eventually undergo a transition to metallic 
and then a superconducting state at pressures consider- 
ably lower than may be necessary for solid hydrogen. Ex- 
ploring the possibility of metallic hydrogen has long been 
a major driving force in high-pressure condensed matter 
science and remains an important challenge in modern 
physics and astrophysics. Recent experimental work on 
SiEU, using diamond-anvil cell techniques, has revealed 
an enhanced reflectivity with increasing pressure [2|, l3[. 
It was found [3[ that solid SiH4 becomes opaque at 27- 
30 GPa and exhibits Drude-like behavior at around 60 
GPa, signalling the onset of pressure-induced metalliza- 
tion. Structural information is the primary step toward 
understanding these observed electronic properties. 

SiH4 has a rich phase diagram with at least seven 
known phases 0, 0, Only one solid phase has been 
reported in the pressure range between 10 and 25 GPa 
and at room temperature, with a monoclinic structure 
(Space group P2\/c) p. The very low hydrogen scat- 
tering cross section of hydrogen-containing materials in 
all diffraction methods makes structural determination 
very difficult, specifically in determining the H positions. 
Although the neutron diffraction is powerful in detecting 
the H-bonding structures, the current accessible pressure 
range of this technique is limited to 30 GPa 0. There- 
fore it is still impossible to use neutron diffraction to 
obtain the interesting structural information of SiEU in 
the metallic state. The challenge of experimentally or 
theoretically determining the high-pressure structures of 
SiHU is still enormous. 

The sequence of SiH4 structures provides the basis for 
understanding whether the material is a favorable candi- 
date of a high temperature superconductor. For a metal- 
lic Pman SiH4 phase, Feng et al. [8] obtained a super- 
conducting transition temperature T c of 166 K at 202 
GPa by using the electron-phonon coupling strength for 
lead under ambient pressure as both materials have the 



same characteristic density of states per volume at rel- 
evant pressures. Pickard and Needs [9[ also studied the 
structural properties of SiEU and mentioned the possibil- 
ity of superconductivity in a C2/c phase. All previous 
work [6|, la, [9[ was done without including the calculations 
of the phonon spectra and electron-phonon coupling pa- 
rameters. Later phonon calculations by Yao et al. [10( 
showed that the Pman structure is in reality not sta- 
ble and that a new C2/c structure is dynamically stable 
from 65 to 150 GPa. This C2/c SiEU phase was pre- 
dicted to exhibit superconductivity close to 50 K at 125 
GPa [TO]. It is not known whether there exists a com- 
mon structural feature that favors superconductivity in 
metallic SiEU, and no study on superconductivity with 
other stable structures has been attempted. 

In this Letter we report a theoretical study of super- 
conductivity in compressed solid SiELi including struc- 
tural, electronic, and vibrational calculations. We find 
six energetically favorable structures in which the PI, 
Cmca, and C2/c structures have layered networks with 
eight-fold SiHg coordination. This layered feature favors 
metallization and superconductivity. The layered phases 
are predicted to have T c 's in the range of 20 and 80 K, 
suggesting that SiEU is indeed a good candidate for high- 
temperature superconductivity. 

To study the structural and electronic behavior of SiEU 
over a wide range of pressure, we used the Perdew-Burke- 
Ernzerhof generalized gradient approximation (GGA) 
density functional and projector augmented wave method 
as implemented in the Vienna ah initio simulation pack- 
age (VASP) [12]. An energy cutoff of 450 eV is used 
for the plane wave basis sets, and 16x16x16 and 8x8x8 
Monkhorst-Pack /c-point grids are used for Brillouin zone 
sampling of two SiH4 molecular cells and four SiH4 molec- 
ular cells, respectively. The lattice dynamical and super- 
conducting properties were calculated by the Quantum- 
Espresso package [l3| using Vanderbilt-type ultrasoft po- 
tentials with a cut-off energy of 25 Ry and 200 Ry for 
the wave functions and the charge density, respectively. 
12x12x12 Monkhorst-Pack /c-point grids with Gaussian 
smearing of 0.02 Ry were used for the phonon calcula- 
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FIG. 1: (color online) The enthalpy versus pressure for com- 
petitive structures of SiH4. The enthalpy of the I4i/a phase 
is taken as the reference point. 
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FIG. 2: (color online) The energetically most favorable struc- 
tures computed for SiH4 polymorphs at various pressures. 



tions at 4x4x4 g-point mesh and double fc-point grids 
were used for calculation of the electron-phonon interac- 
tion matrix element. 

We performed a systematic study of the phase stabil- 
ity of SiEU based on ah initio first-principles calculations. 
Out of more than one hundred structures we studied, six 
new polymorphs of SiH 4 with low enthalpies are found in 
the pressure range from to 300 GPa. In Fig. 1 we plot 
the pressure dependence of their enthalpies along with 
the results for the C2/c and IAi/a structure reported pre- 
viously [9|. A monoclinic structure with P2\/c symmetry 
has the lowest enthalpy below 27 GPa, in good agreement 
with the recent experiments [6]. The P2±/c structure 
consists of four isolated covalently bonded SiHU tetrahe- 
dra with the H atom of one molecule pointing away from 
the H atoms of a neighboring molecule. We find that a 
face-centered orthorhombic structure with Fdd2 symme- 
try appears stable from 27 to 60 GPa. These two phases 
are insulating. 

Near 50 GPa, there are three other competitive, low- 
enthalpy structures with the C2/c, I4i/amd, and IAi/a 
symmetry. Our C2/c structure was based on the UI4 ar- 
rangement, which as exhibits a layered structure. How- 
ever, the C2/c structure predicted by Pickard and Needs 
[9[ forms three-dimensional networks at high pressures. 



For clarification, we name their structure as C2/c(3D) 
and our layered structure as C2/c(2D). There are sub- 
tle differences in band structures between the C2/c pre- 
dicted by Yao et al [10] and our C2/c(2D). Our 
C2/c(2D) structure is composed of six-fold coordinated 
SiH6 octahedra inside a layer, but Si-H bonding between 
different layers is absent. It is still dynamically stable up 
to 250 GPa. However, the stability of the C2/c phase 
of Yao et al [10[ is only stable up to 150 GPa. On 
compression, the six-fold coordinated SiH 6 octahedra are 
transformed into eight-fold coordinated SiHg dodecahe- 
dra with a S2H2 bridge-like bonding arrangement within 
the layer. Above 110 GPa, the C2/c(2D) structure be- 
comes metallic. The IAi/amd SiH4 phase is composed 
of eight-fold coordinated SiHg dodecahedra and every Si 
atom shares two H atoms with other Si atoms, as in the 
Mi /a structure. The Mi/amd-type S1R4 is a semimetal 
between 40 and 70 GPa. 

The I4i /a structure was found to have the lowest en- 
thalpy over a wide pressure range (60 to 220 GPa), consis- 
tent with previous calculations [9] . In an early powder x- 
ray diffraction study [lij], I4i/a was considered as one of 
most plausible structures for the low-temperature phase 
II. Raman measurements [3] indicate that this structure 
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FIG. 3: (color online) (a) Electronic band structure for the 
Cmca SiH4 phase at 250 GPa. The horizontal line shows the 
location of the Fermi level, (b) The phonon dispersion and 
phonon density of states (PDOS) projected on Si and H atoms 
for Cmca SiH 4 at 250 GPa. 



may not exist for SiHU under high pressure and at room 
temperature. SiHU in the IA\/a structure is believed 
to be stable only at low temperature. Between 220-270 
GPa, a metallic Cmca structure has the lowest enthalpy. 
Upon further compression, the C2/c(3D) phase possesses 
the lowest enthalpy between 270-300 GPa, which con- 
firms the previous calculations |9j]. The Cmca phase is 
also a layered structure that consists of eight-fold coor- 
dinated SiHg dodecahedra. In each layer, the Si-H bond- 
ing arrangement is the same as that in the IAi/a struc- 
ture. Thus, the Cmca phase can be viewed as the two- 
dimensional analogue of the IA\/a phase. 

In the 60 to 100 GPa range, we found a metallic 
triclinic structure with PI space group with eight-fold 
SiHg coordination. It has almost the same enthalpy as 
the Cmca structure between 60 and 100 GPa. Both 
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FIG. 4: (color online) Electron-phonon spectral function 
a 2 F{u) vs frequency uj of metallic SiHU with the Cmca, 
C2/c(2D) and C2/c(3D) structures at various pressures. 



of their enthalpies are within 0.25 eV of the insulating 
Mi /a structure. As pressure is increased, the PI struc- 
ture transforms gradually into the Cmca structure. It 
is instructive to note that the PI and Cmca phases are 
metallic over the pressure regime between 60 GPa and 
270 GPa, which is in a good agreement with recent ex- 
periments [3]. Although the insulating IAi/a phase has 
the lowest enthalpy between 60 and 220 GPa, both PI 
and Cmca are good candidates for metallic phases in 
this regime. We thus obtain six energetically favorable 
structures for SiHU at high pressures. The atomic ar- 
rangements for each of these structures is shown in Fig. 
2. The PI, Cmca, and C2/c(2D) phases have layered 
structures and are metallic. 

Figure 3(a) shows the calculated band structure for 
Cmca SiHU at 250 GPa. It can be seen that the Cmca 
structure is metallic. The valence bands cross the Fermi 
level Ep along the YT direction, while the conduction 
bands cross Ep near the Y point. Upon compression, the 
conduction band crossing oss Ep at the V point shifts 
lower in energy, while the valence band across Ep along 
the YT direction only moves up slightly in energy. The 
net effect of the pressure-induced band shifts is to in- 
crease the volume of the Fermi surface and the phase 
space for the electron-phonon interaction. 

The structural stability of each SiH4 phase has been 
examined through lattice dynamics calculations. The 
typical results of the phonon dispersion and projected 
phonon density of states for the Cmca SiHU at 250 GPa 
are displayed in Fig. 3(b). The Cmca stability is con- 
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FIG. 5: (color online) Calculated (a) logarithmic average 
phonon frequency u;i og , (b) electron-phonon coupling param- 
eter A, and (c) superconducting transition temperature T c of 
SiH4 with the Cmca, C2/c(2D), and C2/c(3D) structure as 
a function of pressure up to 300 GPa. 

firmed by the absence of imaginary frequency modes. 
There are weak interactions between the Si framework 
and H atoms over the whole frequency range. The heavy 
Si atoms dominate the low- frequency vibrations, and 
the light H atoms contribute significantly to the high- 
frequency modes. Three separate regions of bands can 
be recognized. The modes for the frequencies below 750 
cm -1 are mainly due to the motions of Si. The bands 
around 200 cm -1 are caused by acoustic phonons. The 
Si-H-Si bending vibrations dominate the intermediate- 
frequency region between 750 and 1200 cm -1 . At high 
frequencies above 1200 cm -1 , the phonon spectrum be- 
longs to the Si-H bond stretching vibrations. 

The electron-phonon spectral function a 2 F(uo) is es- 
sential in determining the electron-phonon coupling 
strength A and logarithmic average phonon frequency 
uo\ og . We have calculated a 2 F(uo) for the metallic SiHU 
phases in the pressure range of interest. Figure 4 shows 
the results for the Cmca, C2/c(2D) and C2/c(3D) struc- 
tures over the pressure range from 70 to 300 GPa. Be- 
low 600 cm -1 the major contributions to a 2 F(uo) come 
from the phonon modes involving Si- Si vibrations, and 
the remaining part of the electron-phonon coupling is 
mainly due to the phonon modes involving H-H vibra- 
tions. The a 2 F(co) on the wide high-energy side is signif- 



icantly higher than that on the narrow low-energy side. 
Thus the high-energy H-H vibrations dominate the total 
A value. Among these metallic structures, the C2/c(2D) 
phase at 250 GPa has a relatively large a 2 F(co) over the 
entire frequency range studied, resulting in a large A. 

We now examine whether superconductivity in metal- 
lic SiH4 is possible, using the T c equation derived by 
Allen-Dynes [15[. In the calculations, we took the 
Coulomb pseudopotential /i* to be 0.1 which was found to 
reproduce T c in MgE>2 [lfj. Figure 5 shows the pressure 
dependence of c^i og , A, and T c for SiH4 with the metal- 
lic Cmca, C2/c(2D) and C2/c(3D) structures. The cal- 
culated o;i g decreases [increases] with pressure for the 
C2/c(2D) [C2/c(3D)] structure. However, A changes 
with pressure in an opposite sense to uo\ og for both struc- 
tures. For the Cmca structure, both co\ og and A do not 
show a monotonic pressure dependence. In these metallic 
phases studied, the variation of T c with pressure is found 
to resemble the A behavior. It is therefore indicated that 
the pressure effect on T c in SiH4 is primarily controlled 
by A. The Cmca phase has a T c of 75 K even at 70 GPa. 
For the C2/c(2D) phase, we calculate a relatively large 
T c as high as 80 K at 250 GPa. A decreasing T c from 47.7 
K at 200 GPa to 26.1 K_at 300 GPa is obtained for the 
C2/c(3D) phase. The PI structure is also estimated to 
have a T c of 46.6 K at 70 GPa. The current results thus 
suggest new possibilities for exploring high temperature 
superconductivity in this hydrogen-rich system. 

In summary, we have investigated the structural sta- 
bility of silane under pressure up to 300 GPa. The P2\/c 
phase is confirmed to be a good candidate for the low- 
pressure insulating phase. Between 27 and 60 GPa, Fdd2 
is predicted to be the structure of another insulating 
phase. At higher pressures, silane enters the metallic 
state having a structure with PI symmetry. As pressure 
is further increased, the PI structure transforms grad- 
ually into the Cmca structure. The three-dimensional 
C2/c structure is most stable only after 270 GPa. The 
layered feature of this material favors metallization and 
superconductivity. The relatively high transition tem- 
peratures in metallic silane are mainly attributed to 
the strong electron-phonon coupling due to the phonon 
modes involving H-H vibrations. 
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